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ABSTRACT 

The inverse cascade of magnetic energy occurs when helicity or rotational instability 
exists in the magnetohydrodynamic (MHD) system. This well known phenomenon has 
been considered as a basis for the large scale magnetic field in universe. However 
nonhelical magnetic energy in a decaying MHD system also migrates toward the large 
scale, which holds vital clues to the origin of large scale magnetic field in a quiescent 
astrophysical system. Zeldovich’s rope dynamo model is considered as a basic and 
symbolistic model of magnetic field amplification. However, the rope model assuming 
specific external forces like buoyancy or Coriolis force is not appropriate for a decaying 
turbulent system without any external force. So we suggest a new dynamo model based 
on magnetic induction equation excluding a forcing source. This model shows the 
expansion and growth of magnetic field (flux) is basically the redistribution of energy in 
the system. The transfer of magnetic energy is in fact a successive induction of magnetic 
field resulted from the interaction between the fluid motion and seed magnetic field. We 
also discuss about an analytic theorem based on the scaling invariant MHD equation. 

Subject headings: galaxies: magnetic fields, decaying turbulence, dynamo 


Introduction 


Magnetic field is a ubiquitous phenomenon in universe which is full of conducting fluids (plas¬ 
mas). The interaction between magnetic field and plasma is known to play a crucial role in the 
evolution of celestial objects like pulsars, jets, galaxy, or GRBs etc. So better comprehension of 
the rudimentary interaction causing the MHD phenomena will more clearly explain the past and 
present universe, and how the universe will evolve. However since the interaction between the 
magnetic field and conducting fluid is basically a nonlinear process, the intuitive understanding 
of MHD phenomena without its exact solution is not easy nor reliable. At the moment the gen- 
eral solution of MHD equation is not yet known. Only a few approximate (stochastic) solutions 
([Kraichnan and Naearaian 1119671 : iKazantsev Ill968l : IPouauet et al Ill976l : iKrause and Radler Ill980lf 
are available in a limited way, which means our comprehension is also restricted as much. So com¬ 
putational simulation is widely used as an alternative method. However the exact interpretation 
of simulation, which is the numerical calculation of MHD equation, requires precise understanding 
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of the analytic theories on MHD. So the investigation of a theory-based intuitive model with the 
numerical simulation is a practical and helpful way to get the point of MHD phenomena. 


We briefly introduce couple of models on the origin of magnetic field in the early universe. We show 
simulation results for the amplification of (nonhelical) magnetic field in a decaying MHD system . 
Then we discuss about a theory derived from the scaling invariant MHD equation (jOlesen 1119971 1. 
Finally we will introduce an intuitive dynamo model for the general mechanism of magnetic energy 
transfer in plasma. They give us clues to the origin of magnetic field in the quiescent astrophysical 
system. 


For the origin of magnetic field roughly two hypotheses are accepted: primordial (top-down) and 
astrophysical (bottom-up) model. Primordial hypothesis supposes magnetic field could be gener¬ 
ated as the confo rmal invariance of electrom agnetic field was broken during the inflationary period 
of early universe (jTurner and Widrow Ill988l b After the expansion, magnetic field could be succes¬ 
sively generated through the cosmic phase transition such as electroweak p hase transition (EWPT) 


or quantum chromodynam ic transition (QCD) from quark to hardron (jGrasso and Rubinstein 


20011 : iSnbramanian 120151 1 . The magnitude of generated magnetic field are thought to have been 
~ 10“®^G on a 1 Mpc comoving scale (durin g inflation), ~ 10“^®G (EWPT), and ~ lO^^'^G 


Bo 

(QGD) on a 10 Mpc scale (jSigl and Olinto lllOOTI ). The correlation lengths of these seed fields were 


limited by the scale of particle horizon 1cm (EWPT) and ~ 10^ cm (QGD). On the other 

hand astrophysical hypothesis, posterior to the primordial inflation, suggests that the seed field 
was generated by the plasma effect in the primeval astrophysical objects like galaxies or clusters. 
The strengths of magnetic seed fields in these models are inferred to be in the range between 
i3o~ 10“^^G (Biermann battery effect) and ~ 10“^®G (Harrison effect). 


However, whether the seed magnetic fields were originated from primordial or astrophysical model, 
the inferred strengths are too weak for the currently observed magnetic fields (~ /uG). Also the 
inferred correlation length, which should be limited to the particle horizon at that time, is too tiny 
compared with the t hat of prese nt magnetic field in space. Definitely the initial seed fields must 
have been amplified (|Gho Il2014l b This process, called dynamo, is essentially the redistribution of 
energy in the MHD system. If the magnetic energy is cascaded toward lar g e scale with its inten¬ 


sity in creasing, this is called large scale dynamo ILSD. iBrandenburg (j200lh . iBlackman and Field 


(j2002l H. In co ntrast if the energy is ca s caded toward small scale, the process is called small scale 
dynamo (SSD. iKulsrud and Anderson (119921 )). Also if the growth rate of magnetic field depends 
on the magnetic resistivity, the process is called slow dynamo; otherwise, it is called fast dynamo. 
The migration and increase of magnetic field are influenced by many factors whose critical condi¬ 
tions are still in debate. But, we will not discuss about the criteria at this time. We focus on the 
qualitative mechanism of energy transfer in the MHD system. 
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As mentioned dynamo is essentially energy redistribution. In the turbulent system the direction of 
energy transfer is related to the conservation or minimization of system variables. For example in 
(quasi) two dimensional ideal hydrodynamics the kinetic energy v? and enstrophy = (V x u)^ 
are conserved. Then using ~ ^ Ev(k,t)dk = 0 and 4a;^ ^ ^ k‘^Ev{k,t)dk = 0 we 
can trace the migration of kinetic energy like below ( Davidson Il2004l i: 

■ kEv{k,t)dk 


^ ^ d / f^kEv{k,t)dk \ 

dv' dt\ ]^ Ev{k,t)dk ) 


A 

dt Jo 


( 1 ) 


f^Ev{k,t)dk 

With the definition of initial centroid of (Ey) hy kc = f kE{k)dt/ f E{k)dt, this result can be 
rewritten like 

1 iJ^{k-kcfEv{k,t)dk 




( 2 ) 


f^Ev{k,t)dk 

Since the naturally spreading energy spectrum Ey{k,t) makes f^(k — kc)‘^Ey{k,t)dk increase, 
the peak of Ey retreats toward the large scale with time; i.e., the inverse cascade of Ey. In three 
dimensional case where the enstrophy is not conserved, Ey migrates toward the smaller scale which 
has larger damping effect. But with the kinetic helicity Hy (= (u • w), w = V x u ~ u), which is 
a conserved variable in the ideal hy drodynamic system, the energy can be inversely cascaded with 
the triad interaction among eddies (jBiferale et al II2012I '). 


In the MHD system without helicity, the kinetic and magnetic energy cascade toward small 
scale. But if the MHD system has a nontrivial kinetic helicity Hy, or magnetic helicity Hm (= 
(A-Bb B = Vx A A), magnetic energy is cascaded inversely (iBrandenburg and Subramanian I 
200a). This phenomenon can be derived using the conservation of energy and helicity, which 


is also related to the minimization of energy in the system. Since the MHD system converges 
into the absolute equilibrium state, the system variables can be generally described by Gibbs 
functional composed of ideal invariants: (tt^ + (A • B), and (u • B). The ensemble average of 
a N-dimensional system variable is represented by the partition function ‘Z’ and the generalized 
coordinate ‘p\ (A) = f pA{p, q) dp^^ dq^^, p = Z~^exp[—W], Z = j exp[—W]dp^^dq^^, 

W = \iE + A 2 A • B + A 3 U • B. However, the system invariant itself does not always guarantee the 
inverse cascade. What determine the direction of cascade are the structures of ideal invariants and 


Lagrangian multipliers: Ai, A 2 , and A 3 (iFrisch et al. Ill975l ). From a fluid point of view, in addition 


to the statistical concept, we can explain the inverse ca scade of Em using the mean field theory 

2ni2al lbll. The profiles of large s cale inagne tic 


(jBlackman and Field 120021 : IPark and Blackman 


energy Em and helicity (A • B) in Fourier space (k = 1) are represented like below (jPark II2014I ): 

d 


dt 


Em = a{A-B)— 2{P+ r])EM, 


m 


(A-B) = 4aEM-2il3 + p){A-B), 


( 3 ) 

( 4 ) 


where 


a = 


((j • b) - (u • u))dT, 13 = ^ J {u^) dr. 


( 5 ) 
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Here, (j-b), (u-a;), and {u^) are respectively small scale current helicity, kinetic helicity, and kinetic 
energy. These equations indicate the energy in small scale is transferred to the large scale through 
the small scale helical field (a (j.b)-(u •a;)), and the dissipation of large scale magnetic field 
is related to magnetic resistivity and small scale kinetic energy (/3 ~ This equation is valid 

for both a driven and decaying helical MHD system. 


Also the cross field diffusive flow (jRvu and Yu Ill998l l or Bohm diffusion effect ([Lee and Rvu 11200711 
was forced to increase the flux of magnetic field. Their concepts are somewhat different from the 
typical dynamo theory, but include the fundamental question on the dynamo. 


In a driven MHD system the saturated field profile is decided by the injection scale, external driving 
force, and the intrinsic properties like viscosity & magnetic re sistivity, not by th e initial conditions 
whose effect disappears within a few simulation time steps ( Park I 2013 . 2014 1. In contrast in a 
free decaying MHD system the profile of decaying field is determined by the initial conditions and 
selective decaying speed besides other intrinsic properties. Especially the selective decaying speed 
due to the different spatial derivative order in the ‘system invariants’ arouses academic interest 
besides many applications: the life tii ne of star-formin g clouds, measurement of galactic magnetic 
field, or fast magnetic reconnection (jBiskamp l2008l ). However, we should note the free decay 
does not mean the interaction between magnetic field and plasma is forbidden. While the total 
energy decays with time, the magnetic field still interacts with the fluid to migrate among eddies. 
Whether or not the MHD system is driven, the helical Em is inversely cascaded as long as ‘a’ in 
Eq.(l3|) is predominant over the dissipation effect. Moreover the magnetic helicity in large scale, the 
statistical correlation between different compo nents of magnetic fields, is resident to the turbulent 
diffusion so that the decay speed slows down ([Blackman and Subramanian Il2013l b On the other 
hand the inverse transfer of Em i n a decaying system without he l icity or shea r is contradictory to 


the t ypical MHD dynamo theory (jOlesen 


1997: Shiromizu 


1998; 


Zrake II2014I : (Brandenburg et al 


20141 ). This phenomenon implies the energy released from the past events like supernovae could be 


a source of large scale magnetic field structure in a quiescent astrophysical object which does not 
have a significant driving source. 


Usually the amplification of magnetic field (fl ux) due to the fluid mo tion is explained by Zeldovich’s 
rope model ‘stretch, twist, fold, and merge’ (IZeldovich et al Ill983l l. However, since the model as¬ 
sumes the influence of external force, it is not appropriate to the growth of large scale magnetic 
field in a decaying turbulence. 
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1. Simulation and method 


We used PENCIL CODE (jBrandenburg Il200lh for the weakly compressible fluid in a periodic 
box (Stt^). MHD equations are basically coupled partial differential equations composed of density 
‘p’, velocity ‘u’, and vector potential ‘A’ (or magnetic field B = V x A). 


-pV ■ u, (6) 

—Vln /?+-(V xB) xB + i/(V^u + - VV • u), (7) 

p 3 

uxB —T/VxB + f. ( 8 ) 

Here D/Dt{= d/dt + u ■ 'V) is Lagrangian time derivative to be calculated following the trajectory 
of fluid motion. In simulation we used 0.015 and 2 x 10~® for the kinematic viscosity V’ and 
magnetic diffusivity ‘t/’ respectively. This setting is to realize the large magnetic prandtl number 
Pf'M ^=v/r]) in the early universe as similar as we can. The function f(x,t) generates the random 
nonhelical force with the dimension of electromotive force (EMF, u x B). The nonhelical magnetic 
energy maximizes Lorentz force J x B so that the fluid motion in plasma is efficiently excited. 
Turning on and off the function (0 < t < 1, simulation time unit), we imitate a celestial system 
which was left decaying after it had been driven magnetically by an event in the past. 


2l 

Dt 

Du 

15t 

dA 


2. Results 

Fig jl(a)| shows the evolution of kinetic and magnetic energy spectrum in a decaying MHD 
system after the initial nonhelical magnetic forcing (0 < t < 1). Black dotted line indicates kinetic 
energy spectrum Ey, and red solid one indicates magnetic energy spectrum Em- The increasing 
thickness of line means the lapse of time from t = 5 (thinnest) to t = 1342 (thickest). The overall 
Ey and Em decrease after the forcing stops, but each spectrum does not monotonically decay. For 
5 < t < 60 Ey in the small scale {k >~ 25) grows before it turns into the decay mode at t ~ 60, 
and the lagging decrease of Ey in the large scale gives some hints of the influence of eddy turnover 
time ‘r’. The profile of Em shows the similar, but somewhat different pattern. Em in small scale 
grows faster than that of large scale, surpasses its initial magnitude, and then begins to decrease 
at t ~ 60. In contrast, Em in large scale keeps growing far longer before it eventually decays. The 
overall evolution of magnetic field lags behind that of velocity field, which is a typical feature of a 
large Pvm MHD system 


Fig jl(b)| in fact has the same information as that of Fig jl(a)[ But it shows the important features 
of a decaying nonhelical field more clearly. After the event stops at t ~ 1, both Ey and Em grow 
temporarily before they decay. The onsets of large scale Em and Ey lag behind those of small 
scale energy but keep growing longer, which implies the influence of r, Pvm, and energy supply on 
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the evolution of large scale field. The advection term u • Vu’ or pressure Vp’ transfers Ey to 
make the system homogeneous and isotropic without the help of Em- But for the transfer of Em, 
Ey is a prerequisite as the source term ‘V x (u x b)’ in the magnetic induction equation shows. 
The energy in the kinetic eddy is non locally transferred to the magnetic eddy through B • Vu, and 
the magnetic energy is locally transferred through -u • VB. Therefore the long lasting Ey in large 
scale can induce the migration of Em toward the large scale, which also keeps longer due to the 
proportionally increasing eddy turnover time (~ 1/A:) and decreasing magnetic diffusivity (~ k"^). 
The inverse transfer of Em without helicity or shear deviates from the principle of typical dynamo 
theory. However the seemingly inconsistent phenomenon is true and implie s the fundamental prin¬ 
ciple of energy migration in the magnetized plasma system. lOlesen (jl997l l suggested the (inverse) 
transfer of energy be an essential phenomenon in the scaling invariant (self-similar) MHD equations, 
and the direction of energy transfer should be decided by the initial distribution of given energy 
^E{k,0) ~ k'^’. The self similarity theorem shows that energy is inversely transferred if (7 < —3; 
otherwise, the energy cascades forward. The inverse transfer of Em in a decaying MHD system im¬ 
plies that an event, which emitted electromagnetic energy in the past, can be an origin of large and 
small scale magnetic field observed in the present universe. As an application of this plot, an ob¬ 
server at ‘t ~ 1000’ can find large and small scale Ey & Em without any driving source nor helicity. 


Fig,2(a) 2(b) show the decaying Ey and (helical) Em in Fourier space. Fig 2(c) includes the 


evolving profiles of total Em solid red line), helical Em {{k A ■ B)/2, dotted red line), and 

Ey (dashed black line) in real space. Fig j2(d) shows the evolving helicity ratio EHm/^-Em' at k =1, 
5, and 8 . The conditions are the same as those of nonhelical case except that the system is initially 
driven by the helical magnetic energy. The helicity makes some distinct features discriminated 
from the nonhelical decaying turbulence. The migration of “'Em peak’ appears clearly, the strength 
of kinetic energy and overall decay rate of energy are lower than those of the nonhelical energy 
system. Like the nonhelical MHD system the magnetic energy is transferred to the kinetic eddy 
through Lorentz force. In principle the fully helical magnetic energy cannot be transferred to the 
kinetic eddies. However the unstable turbulent motion lacking in memory effect generates nonhelical 
component so that kinetic eddies receive energy from magnetic eddies. So although the strength 
of Ey is weaker, there is limited contribution of Ey to the induction of nonhelical magnetic field. 


Careful comparison of Fig 1(b) 2(c) 2(d) shows that the kinetic energy is transferred toward both 


directions; consequently, Ey in large and small scale leads to the forward and backward migration 
of Em- In the early time regime (1 < t <~ 50), the helicity ratio of magnetic energy keeps 
relatively constant in spite of the elevation of Em- When the helicity ratio in large scale Em begins 
to accelerate {t >~ 50, Fig 2(d)), that of small scale Em starts falling. We can infer the transfer 
of total Em due to Ey precedes that of helical Em due to the small scale helical magnetic field 
{‘a effect). And then ‘a coefficient’ interacts with large scale magnetic held directly leading to the 
cascade of small scale helical magnetic energy inversely. The role of pressure and advection term 
seems to be independent of the helicity ratio in the system. 
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3. Theory: Models of energy transfer 
3.1. Inverse cascade of helical magnetic energy 


As Fig 1(a) 2(d) show, the magnetic energy in a decaying MHD system can be transferred 
toward the larger scale whether the field is helical or not. In case of helical field, toroidal mag¬ 
netic field interacts with the helical fulid motion to generate (amplify) poloidal magnetic held, and 
the poloidal component generates (amplihes) the helical component through the interaction with 
plasma. The effect of helicity in small scale is represented by ‘o’, which can be considered as 
an independent coefficient in a homogeneous and isotropic (without rehection symmetry) system. 
According to the relative contribution to the held amplihcation, the helical dynamo is divided into 
dynamo’, ‘afl dynamo’, dynamo’ (Here ‘fl’ indicates the differential rotation effect). 


The large scale magnet ic energy an d helicity in a decaying system are represented by the solutions 
of coupled Eq.([3]), ([4]) ( Park 2014 ): 


2HMit) = (Fm(0) + + (//^(o) - 2F;M(0))e-2/o(“+/5+^)'^-, (9) 

4EMit) = (Hm(0) + 2EMiO))e‘^^oi»-^-v)dr _ (h^(o) _ 2EMiO))e-^^o('^+^+v)dr_ (^ q ) 


(Here Em{0) and Hm{0) are the initial large scale magnetic energy and helicity.) 

With the positive initial magnetic helicity (a > 0), the hrst terms on the right hand side (RHS) 
are dominant. However, since there is no driving source, ‘a ( (j • b) - (u •w))’ and ‘/3 (~ {n^)y 
eventually decay and converge to ''zero'. But the combined index 'a — /3 — rj' decreases to become 
negative with the hnite diffusivity ‘ 77 ’. Consequently the strengths of (large scale) magnetic energy 
and helicity increase hrst, reach the peak, decay and converge to zero eventually, which is consistent 
with the held evolution shown in Fig )2(c)[ If the initial magnetic helicity in small scale is negative, 
i.e., a < 0 the second terms in RHS decide the prohles of EM{t) and HM{t) with the negative sign. 


3.2. Inverse transfer of nonhelical magnetic energy 


The inverse transfer of decaying nonhelical magnetic energy cannot be explained by Eq. ©, dini) 
because of the negligible ‘o’. Instead there have been trials to explain the pheno menon qualita¬ 


tively using the scalin g invariant Navier Stokes a nd magnetic induction eq uation (jOlesen 1119971 : 


Ditlevsen et al. 

2004 

) with numerical simulation ( 

Brandenbure et al 

2014) 


free decaying MHD equations 

9u 

'dt 


= -u • Vu-h b • Vb - VP-Fz/V u, 

'-V-' 

jxb—Vp 

= b • Vu — u • Vb-l-r? V^b. 

Vx (uxb) 


( 11 ) 


5b 


( 12 ) 























are invariant under the scaling transformation: r —)• Zr, t ^t, u —>■ Z^u, u — >■ b Z'^b, 

r] —>• Z^'*“^r/, P —)■ Z^^P, where ‘Z’ and ‘Zi’ are arbitrary parameters. Then the scaled kinetic energy is 


Ey(A:/Z,Z^ ’^t,Ll,K/l) = f d^xd^y ^^(u(Zx,Z^ ^t)u{ly,P ^t)) 

(2ttY J2^/k 




(13) 


Similarly the scaled magnetic energy is represented like ‘EmCZc/Z, LI, K/l) = Z^'''^^Em(Zc, t, L, K)\ 
So the energy density spectrum is ^Ev,M{k/l,l^~^t,lL,K/l) = M(k,t, L, K)\ wh i ch ca n 

be defined as ^Ey luik, t) = k~^~‘^^ijjv^M{k, t)' with an arbitrary function ^Y>{k, Z)’ (jOlesen (jl997l '). 

(120041 1 ■ references therein). These scaled energy density relations lead to 


Ditlevsen et al. 


Yv,M{k/l, P ^Z) = V’V',m(Zc, Z). 


(14) 


If we differentiate Eg. (1141) with respect to ‘Z’ and then set ‘Z = 1’, we derive a differential equation: 


+ ii - h)t?E^ = 0^ 


dk 


dt 


(15) 


The general solution of this equation implies ^Yv,M{k, Z)’ is the function of which implies the 

inverse transfer of energy with increasing time ‘Zc ~ For example if the primordial energy 

is ^Ey^Mik, 0) = (q = —1 — 2Zi)’, the energy at ‘Z’ is ^Ey^M{k, Z) = k~^~‘^^Yv, M{k^~^t) = 

k^'py^M{k^^^^'>^‘^ty. More clearly, when Zi < 1 (i.e., q > —3), the decaying energy can be transferred 
inversely. Briefly the energy migration is the result of scaled energy density relation, which is the 
intrinsic property of MHD equations. However this criterion is not valid in the whole range. The 
integr ation of Ey_ m ^ k^ Yv to get the total energy is known to yield an inconsistent 

result (|Ditlevsen et al. 1120041) . The saturated energy spectrum of Ey and Em in a ‘large Pvm MHD 
syste m’ driven by the nonhelical kinetic energy has the relation of = k'^Ey' in the subviscous 
scale (jPark l2015l l. When the background magnetic field Lext is strong, turbulence due to the 
nonlinear effect becomes relatively weak compared with that of the guiding background m agnetic 
field. This leads to the energy spectrum Ey ~ and Em ~ k~^ (jLazarian et al.l 120041) . If the 
forcing stops, Ey should be cascaded toward small scale but Em should be transferred toward large 
scale. However Ey and Em cannot migrate oppositely. With the strong b^xt the assumption of 
self-similarity in the MHD equations is not valid. Moreover it is not yet clear if the scaled energy 
density relation can be applied to the helical field MHD system. More study is necessary. 


3.3. New model of energy transfer based on MHD eqnation 

We suggest a new dynamo model for the transfer of Em- This model is not limited to the 
free decaying turbulence; rather, it covers the general mechanism of converting kinetic into mag¬ 
netic energy and the local magnetic energy transfer. As mentioned dynamo is in fact the result 
of interaction between the velocity field ‘n’ and magnetic field ‘b’, i.e., EMF (u x b) and cross 
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helicity (u • b). While EMF plays a role of explicit source of th e magnetic field, the cross helicity 
constrains or suppresses the field profile implicitly (jYokoi Il2013l l. In this paper we do not consider 
the effect of cross helicity; instead, we show how the curl of EMF V x (u x b) and its decom¬ 
posed terms ‘b • Vu’, u • Vb’ induce the magnetic field. The MHD system is filled with various 
kinds (scale, magnitude) of magnetic and velocity fields. If both velocity and magnetic field are 
composed of a single field vector, ‘V’ or ‘Vx’ is meaningless and no dynamo occurs. The fields 
should be at least locally inhomogeneous and anisotropic although the macroscopic system may 
be homogeneous and isotropic. So we consider two simplest cases that show the physical process 
of dynamo clearly: ‘u’ (or b) is a single field vector, but ‘b’ (or u) is composed of plural field vectors. 


In Fig, 3(a) we assume the magnetic field heads for the gradient of velocity field. That is, we set 
the magnetic field: ‘b = (0, 0, 6)’ and velocity field: ‘uj = (0, Ui{z), 0)’. The strengths of velocity 
fields ‘Uj {i = 1, 2)’ and ‘U’ are assumed to be the function of ‘z’ in order of U 2 < ui < U. Also 
the velocity field is not too strong compared with ‘6’. Then we can make clear the meaning of first 
term ‘b • Vu’ in Eq. (|12p . Mathematically it is the contraction of second order tensor into the first 
order one, i.e., vector: 


b ■ Vu bz ■ 


^ d ^ d ^ d \ , . ^ , duAz) ^ 

+ 2^0“ + p*(^)y = 0—^— y- 

ax ay oz J oz 


(16) 


Since strength of ‘uj(z)’ increases toward ‘i’, the induced magnetic field ‘bj^^’ is parallel to ‘uj’. 
And ^hind is merged with ‘by’ which is partially dragged ‘b’ by the fluid motion ‘uj’. If we assume 
‘b||’ is the strongest at ‘f7’ and ^duijdz'’ is uniform, the transferred magnetic field ‘6md+&||’ will have 
the largest value at ‘17’. However since ‘b- Vu’ is the part of V x (u x b), some physical information 
is missing. So we need to derive the induction of magnetic field through the curl of EMF directly. 
As mentioned ‘b’ interacts with the fluid ‘uj’ to generate EMF varying from the smallest (u 2 x b) 
to the largest (U x b). These differential magnitudes in EMF yield the clockwise rotational effect 
V X (uj X b) > 0, which induces the magnetic field h^d ~ / drV x (u* x b). Therefore, hind and 
b||_ i are in the complementary relation to cause the growth of net magnetic field. If the strength of 
‘b||’ is assumed to be uniform, bind + b\\^ i will be the strongest at U 2 - However if ‘6||’ is proportional 
to the strength of ‘uj’, the combined magnetic field will have the largest value between ‘uj’ and 
‘U’. This can be one of the reasons why the peak of Em in a non helically driven MHD system 
(small scale dynamo) is located between the injection and viscous scale. Of course the influence of 
eddy turnover time ‘r*’ and dissipation effect k'^ bi need to be considered in order to explain the 
reason more exactly. On the other hand if the direction of b is reversed (b • Vu < 0), hind is in 
the opposite direction of by resulting in the dissipation of net magnetic field. This is the process 
of nonlocal energy transfer from Ey to Em- 


In the same way, we can explain the local energy transfer in magnetic eddies as shown in Fig 3(b)[ 
At this time the fluid ‘u’ is heading for the decreasing egui-magnetic field line (u • Vb < 0). Each 
velocity and magnetic field can be represented like ‘u = (0, 0, u)' and ‘b* = (0, bi{z), 0)’. The 
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strength of magnetic field is the function of ‘z’ in the order of 62 < 6 i < -B and dhi{z)ldz < 0 . 
Then the second term u • Vbj’ in Ea. ()12p is 

/a d d 

-uz ■ \ x-^ + y-^ + z^ 

\ ox ay oz 

So the direction of induced magnetic field ‘'hind is y, which is parallel to ‘bj’. Also their sum¬ 
mation is the largest at B. However, this cannot explain local energy transfer. So we need to 
derive the induced magnetic held using the mathematical dehnition of curl operator again. The 
interaction between ‘u’ and ‘b’ generates EMF varying from the smallest (u x b2) to the largest 
(u X B). Their differential strengths create the rotational effect, which induces magnetic held 
hind ~ f dr's/ X (u X bj). This induced ‘bj„^’ is parallel to ‘bj’ and largest at ‘b2’ indicating the 
energy migration from the strongest magnetic held to the weakest one, which is independent of an 
eddy scale. This explains the inverse transfer of Em in a decaying MHD system. Like the nonlocal 
energy transfer the direction of local Em transfer is determined by the location of velocity held, 
relative strength of energy among eddies, eddy turnover time, dissipation effect. Also we can see 
that ‘-hind generated by the oppositely directed ‘u’ annihilates ‘bj’. 


Oi{^)y = -u y 


u ■ 


dbi{z) 


dz 


(17) 


The induced magnetic held also constrains the huid motion. The current density ‘J’ in Lorentz force 
lEo llip is derived from the how of electric charge ‘gu’ carried by the huids. In view of macroscopic 
huid model, ‘ J’ is replaced by V x B, meaning the induced magnetic held is also able to constrain 
the huid motion. In Fig, 3(a) where the magnetic held ‘bz’ crosses the velocity helds, the directions 
of hind and V x hind are ‘y’ and ‘x’ respectively, the direction of (V x bj^^) x hind is ‘E. This 
appears to perturb the huids ‘u*’ perpendicular. However, (V x hind) x bind can be decomposed into 
the magnetic pressure and magnetic tension Vbj„d’. Magnetic tension cancels the 

magnetic pressure parallel to the held so that x bj^^’ seems to press the huid perpendicular. 
However, this tension is actually hindhind ■ ^b\^j2 -b bf^^hind ■ Vbi„rf = bj^^V \\bind/‘^ + bfnd^ where 
‘k’ is a measure of the curvature of hind- So the exact Lorentz force is —V_l6^^^/2 + bfndf^- This 
curvature related force gives energy to the huid motion. Thermal and magnetic pressure suppress 
the huid motion, but magnetic tension boosts it. On the other hand in Fig 3(b) (V x hind) x hind 
perturbs the huid motion. 


These two cases explain the amplihcation, more exactly migration of Em and its constraint on 
Ey. Real distribution of Ey and Em is more complicated, but basically they can be replaced by 
the combination of these two structures. What we have neglected is the effect of cross helicity 
(u • b). When magnetic held ‘b’ is (anti) parallel to the huid motion ‘u’, EMF is zero so that 
‘hind^ is not generated. But more detailed stochastic analysis shows the effect of cross helicity is 
generated when the fourth order moment is decomposed into the combination of second order one 
to close the MHD equation s : (uu bb jnd) — ( uu)(b bi„d) + {ub){ubind) (Quasi Normal approximation. 


Kraichnan and Naearaian ( 19671 1. IPark I ( 2015l lL which qualitatively matches the cross helicity 


term ‘(uj-bj^^)’ shown in Fig 3(a)|or[3(b)[ For detailed investigation of cross helicity and anisotropy 
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in the magnetized plasma system, more elaborate method is required. But we will not discuss about 
those topics here. 


4. Summary 

We have shown the inverse transfer of Em in a free decaying MHD turbulence. In the typical 
2D hydrodynamic and 3D MHD turbulence system, the inverse transfer of energy is related to 
the conservation of physical invariants such as energy, enstrophy, or magnetic helicity. The overall 
mechanism of inverse cascade of helical Em is well understood. However it is still tentative if the 
magnetic helicity is prerequisite to LSD. Moreover since the helicity is negligible or completely 
absent from some celestial objects that are filled with various scales of magnetic fields, the exact 
comprehension of (inverse) transfer of general Em is important to understand the origin and mech¬ 
anism of magnetic field evolution in the present universe. 


As we have seen, the inverse transfer of nonhelical Em is the result of intrinsic property of scaling 
invariant AIHD equations. The evolving energy spectrum at ‘t’ is represented by the initial energy 
distribution Ey^Mik, 0) ~ k'^ and However this relation, limited to the local range, 

may induce some inconsistent inference for the direction of energy transfer. 


We introduced a dynamo model based on the magnetic induction equation dh/dt ~ V x (u x b) ~ 
b • Vu — u • Vb. This model is not limited to the inverse transfer of decaying Em, but it explains 
the general mechanism of Em transfer. The migration of Em is in fact the continuous induction 
of magnetic field guided by Ey- The energy in the kinetic eddies is non locally transferred to the 
magnetic ones through ‘b • Vu’; and, the energy in the magnetic eddies is locally transferred to 
their adjacent ones through u- Vb’. As the model shows, what decides the direction of magnetic 
energy transfer is not the eddy scale size, but the guide of velocity field and relative energy differ¬ 
ence between eddies. In a mechanically driven system, the pressure and advection term transfer 
Ey chiefly toward the smaller scale whose eddy turnover time decreases by ~ 1/k but dissipation 
effect elevates by ~ A:^. Then Ey cascaded to the small scale induces Em continuously forming 
the peak of Em in small scale. But strictly speaking energy migrates toward both directions. If 
the forcing in the system stops, Ey in small scale fades away more quickly than that of large 
scale. Consequently Ey in large scale generates Em in large scale and decides the profile of energy 
spectrum in the system. Again magnetic back reaction due to Lorentz force generated by this Em 
constrains (suppresses or boosts) the fluid motion, which is consistent with other numerical results. 


In addition to this current model and numerical simulation, a more detailed analytic theory and 
various numerical simulations with arbitrary helicity ratio are required. Since we need to trace 
the evolution of Ey and Em together, Eddy Damped Quasi Normalized Markovianized (EDQNM) 
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approximation is a more suitable meth od than other stochastic models ([Kraichnan and Naearaian 


19671 : IPouauet et al 1119761 : iPark Il2015l h According to the numerical test of EDQNM for a unit Pvm 
MHD system! Son Ill99dh the spectrum of low k is not modified, but the peaks of Ey and Em clearly 
migrates toward large scale. We think the approximation of EDQNM with large Ptm can show 
clearer transfer of Em toward the large scale and give us the hints about the exact representation 
of '>Pv,M{k, t) beyond the inertial range. We leave these topics for the future work. 
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0.1 1.0 10.0 100.0 1000.0 10000.0 
time 


(a) 


(b) 


Fig. 1.— (a) Free decaying Ey and Em spectrum after nonhelical forcing. Vertical axis means 
the energy, and horizontal axis means the wave number, (b) The time evolution of Ey and Em- 
Vertical axis is the energy, and horizontal one is time. The lines of k=l, 5, and 8 represent the 
large, injection, and small scale energy spectrum respectively. The plot shows the increasing energy 
in large and small scale is in fact the transferred one from the injection scale at k = 5. 
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k (t = 5, 60, 350, 700, 1342) 


(a) 



0.1 1.0 10.0 100.0 1000.0 10000.0 
time (helical forcing ^ free decaying) 



(b) 

<kA-B>/<B-B> 



0.1 1.0 10.0 100.0 1000.0 10000.0 
time (helical forcing ^ free decaying) 


(c) 


(d) 


Fig. 2.— (a) Spectrum of decaying kinetic and magnetic energy, (b) Spectrum of decaying total 
magnetic energy {B‘^)/2 and helical magnetic energy {kA. • B)/2. Statistically helicity is the cor¬ 
relation between different components of fields ~ {BiBj). Magnetic energy includes helical 
and nonhelical component, (c) Inverse cascade of total magnetic energy (solid red line), helical 
magnetic energy (dotted red line), and kinetic energy (dashed black line) of fc=l, 5, 8. (d) Helicity 
ratio kHM 

























- 16 - 




u ^ b 


Fbt 

B- 


(•)(ux B) 


> 


bi 

^2 


® (u >^b i) 

(«) (^ X b^) 


y 

11 





!▼ bind~JdTFx (u Xbj) 

y 

u 


(b) 

Fig. 3.— (a) The energy in kinetic eddies is cascaded toward magnetic eddies in a 
(b) The energy in a magnetic eddy is transferred toward its adjacent eddy. 


nonlocal way. 




















